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Molecular dynamics simulations of a short-chain polymer melt between two brush-covered surfaces 
under shear have been performed. The end-grafted polymers which constitute the brush have the 
same chemical properties as the free chains in the melt and provide a soft deformable substrate. 
Polymer chains are described by a coarse-grained bead-spring model with Lennard- Jones interactions 
between the beads and a FENE potential between nearest neighbors along the backbone of the 
chains. The grafting density of the brush layer offers a way of controlling the behavior of the 
surface without altering the molecular interactions. We perform equilibrium and non-equilibrium 
Molecular Dynamics simulations at constant temperature and volume using the Dissipative Particle 
Dynamics thermostat. The equilibrium density profiles and the behavior under shear are studied as 
well as the interdigitation of the melt into the brush, the orientation on different length scales (bond 
vectors, radius of gyration, and end-to-end vector) of free and grafted chains, and velocity profiles. 
The viscosity and slippage at the interface are calculated as functions of grafting density and shear 
velocity. 



I. INTRODUCTION 

Grafted polymer layers have received abiding attention 
because coating surfaces with polymer brushes offers an 
opportunity to control static and dynamic surface proper- 
ties. Polymer brush coatings have been utilized in a wide 
variety of applications ranging from the stabilization of 
colloidal suspensions (e.g., paint), the reversibly tuning 
of wetting and adhesion properties to lubrication, friction 
and wear. 1 For instance, polymer brushes are renown for 
their ability to drastically decrease the friction between 
two solid surfaces. 2 The work of adhesion and wettabil- 
ity are also changed when a solid surface is coated with 
end-grafted polymer layers. Thus, polymer brushes and 
cross-linked rubbers find applications as adhesives. 

Polymer brushes are also interesting from a fundamen- 
tal point of view because many properties are determined 
by the conformational entropy of the extended macro- 
molecules, in particular, the stretching of the polymers 
away from the grafting surface 3 in response to increasing 
the grafting density or swelling the brush with a solvent. 
Increasing the grafting density is an experimentally con- 
venient way of altering the surface properties without 
altering molecular interactions. This gives rise to a rich 
wetting behavior even in apparently simple model sys- 
tems like a polymer drop on top of a brush of identical 
chains. 4-8 If there are few chains grafted onto the surface, 
an increase of the grafting density will result in a smaller 
contact angle because the melt chains benefit from the 
additional attraction provided by the brush. At inter- 
mediate grafting densities this might result in complete 



wetting of the melt on the top of brush. If the grafting 
density is too high, however, the polymer melt will dewet 
from the brush - a phenomenon termed autophobicity. 
Such "entropic dewetting" has been observed experimen- 
tally in polymer networks 9 and polymer brushes. 10,11 Au- 
tophobicity can be traced back to the formation of an in- 
terface between the strongly stretched chains in the poly- 
mer brush and the free chains of the melt. 4 ~ 7 

The interface between brush and melt has many im- 
portant applications: Such interfaces commonly occur 
between two immiscible homopolymers that have been 
reinforced by a diblock copolymer. The diblock adsorbs 
at the interface as to place its blocks into the respec- 
tive phases. Thereby, it lowers the interfacial tension 
and, for sufficiently large surface excess, each block forms 
a polymer brush in contact with the homopolymer-rich 
bulk phase. 12,13 The interdigitation between brush and 
homopolymer chains determines many mechanical prop- 
erties of the interface. 14,15 While the equilibrium prop- 
erties of these copolymer brushes at the interface be- 
tween incompatible homopolymers (as well as brushes 
formed from end-adsorbed chains) resemble those of ir- 
reversibly grafted polymer brushes, their dynamic prop- 
erties might substantially differ due to the in-plane mo- 
bility. Another realization of a brush-melt interface oc- 
curs when a polymer melt dewets from a dense polymer 
brush (or rubber 9 ). These soft deformable substrates 
might offer great control over flow properties. As we shall 
demonstrate below, changing the grafting properties one 
can tune the hydrodynamic boundary conditions. If the 
grafting density is high enough for autophobic dewetting 
to occur,, additional mechanisms of dissipation for mov- 
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ing droplets like lifting of the three-phase contacts line 
(wetting ridge) 16 to balance the normal components of 
the forces on the contact line or chain pull-out 17 can be 
envisaged. 

The equilibrium properties of polymer brushes in var- 
ious solvents have been studied by Monte Carlo 18 ' 19 and 
Molecular Dynamics (MD) simulations 20-23 as a function 
of chain length and grafting density. There are fewer 
studies investigating the equilibrium properties of the 
brush-melt interface. Grest studied the limit in which 
the melt chains are much shorter than the chains of 
the brush. 21 Daoulas et al. performed detailed atom- 
istic Monte Carlo simulations to analyze the equilibrium 
structure a chemically realistic united atom model for 
polyethylene melt in contact with a brush grafted onto a 
carbon substrate. 23 These brush-melt interfaces of long 
entangled chains have also been investigated in the con- 
text of adhesion. Sides and Grest used glassy brush-melt 
interfaces to study chain pull-out and craze formation 
under tensile stress by means of simulations. 24 

Although the molecular mechanisms and the concomi- 
tant boundary conditions (amount of slip) of flow at poly- 
meric interfaces have attracted abiding interest, it is only 
incompletely understood. Many distinct physical pro- 
cesses are involved such as wetting, roughness and bubble 
formation. 25 An accurate prediction of non-equilibrium 
properties is important for the design of nano- and mi- 
crofluidic devices characterized by a large surface-to- 
volumc ratio. In particular, the boundary condition of 
polymer melts between walls has been extensively stud- 
ied experimentally, 26,27 by computer simulation 28,29 and 
by analytic theory. 30,31 The boundary condition and fric- 
tion mechanisms have also been studied experimentally 
with near field laser velocimetry techniques, including 
anchored chains in the solid surface. 27 The sliding of 
two brushes and the non-equilibrium behavior of brush- 
brush interfaces can be studied by surface force appa- 
ratus (SFA) measurements 2 and the experimental results 
have been compared with computer simulations. 21,22,32 ' 33 
These simulations emphasize the correlation between the 
non-equilibrium shear properties and the interdigitation 
and orientation of the two apposed brushes. 

In this work we study the behavior of the brush- 
melt interface of short non-entangled chains of a coarse- 
grained model in equilibrium and under shear as a func- 
tion of the two experimentally accessible parameters: 
grafting density and shear velocity. The grafting den- 
sity is systematically varied over a wide range from the 
mushroom-regime to high grafting densities that corre- 
spond to autophobic behavior. MD simulations in junc- 
tion with the Dissipativc Particle Dynamics thermostat 
are employed to account for hydrodynamic correlations 
without including an explicit solvent. Unlike most calcu- 
lations in the literature we study polymer chains in a poor 
solvent. The attraction between monomers allows us to 
extend this model to studying liquid- vapor interfaces and 
the wetting of a melt on top of a brush of identical chains. 
The remainder of our manuscript is organized as follows: 



In section II we introduce our coarse-grained model and 
the simulation technique. In section III the results of the 
simulations for the structure, orientations and boundary 
conditions of this melt-brush interface in equilibrium and 
under shear are detailed. The last sections contains a dis- 
cussion of our results and an outlook on future work. 



II. MODEL AND THERMOSTAT 

We consider a widely utilized, coarse-grained bead- 
spring model for polymers. 34 This model has been used 
for a variety of thermodynamic conditions, chain lengths 
and physical regimes such as glasses, melts and dilute 
solutions. 35 The potential between neighboring beads 
along the same polymer is modeled by a finite extensible 
non-linear elastic potential (FENE) , 
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where Rq = 1.5a, the spring constant is k = 30e/cr 2 , and 
r = |rj — rj\ denotes the distance between neighboring 
monomers. Excluded volume interactions at short dis- 
tances and van-der-Waals attractions between segments 
arc described by a truncated and shifted Lcnnard-Jones 
(LJ) potential: 



U(r) = U hJ (r)-U hJ (r c ), 



with 
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where the LJ parameters, e = 1 and a — 1, define the 
units of energy and length, respectively. Ui,j{r c ) is the LJ 
potential evaluated at the cut-off radius of twice the min- 
imum of the potential: r c = 2 x 2a. Often a completely 
repulsive LJ potential, r c = 2 e , has been utilized for com- 
putational efficiency. 32 ' 33,36 This choice corresponds to a 
polymer brush in a good solvent. Here, we include attrac- 
tive interactions between segments and, thus, are able 
to describe liquid-vapor coexistence, wetting and droplet 
formation below the 0-temperature, fcgS/e ~ 3. 3. 37 

The substrates at z = and z = D are modeled as 
idealized flat and impenetrable walls which interact with 
the polymer segments via an integrated Lennard-Jones 
potential, of the form 



V^{z) = \A\(^f 

where A = 3.2e is sufficient to make the liquid wet the 
bare substrate. 37 

The lateral positions of the grafted chain heads are 
randomly chosen with a total number of chains given by 
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the desired grafting density. The distance of the first 
bead from the grafting surface is Z\ = 1.2c. 

In our model wall roughness is completely disregarded. 
The friction between the melt of free chains and the sur- 
face stems from the grafted chains only. Our goal is to 
analyze the experimentally controllable properties of the 
substrate provided by the grafted polymer layer without 
any additional source of roughness arising from a possible 
corrugation of the substrate. 

We used a Dissipative Particle Dynamics thermostat 
(DPD) to simulate at constant temperature 38 ' 39 and ac- 
count for hydrodynamic interactions due to the conser- 
vation of total momentum. The equations of motion are 
formally the same as used in Langevin dynamics: 40 ' 41 



Pi = Fj + Fj + Fj (5) 

where ¥i D and Fi R denote the dissipative and random 
forces respectively, rm is the mass of the monomers which 
is set to unity for all the simulations. The difference be- 
tween DPD and Langevin thermostats is that the random 
and dissipative forces are applied in a pair- wise form, 
where the total forces acting on a pair are equal to zero 
and, thus, momentum conservation is obeyed. The ex- 
pression for the forces are the following: 

F >° = E F £ ! F £ = -7^(r«)(^ ' v«)*« (6) 
V." X F£ ; Fg=C W R (r«)¥i;. (?) 

where for each vector a we use a notation a^- = a, — a,, 7 
is the friction constant and ( denotes the noise strength. 
Friction and noise, 7 and £, obey the relation ( 2 = 2IcbTj 
and the concomitant weight functions satisfy fluctuation- 
dissipation theorem, [ui R ] 2 = uj d . 9ij is a random vari- 
able with zero mean and second moment (6ij(t)6ki(t )) = 
(SijSji + SuSjk)S(t — t ) and the weight functions are cho- 
sen to be: 

[^] 2 =^={ (1 - r/rc)2 ' r<rc ' (8) 
L J \0, r > r c 1 ; 

which is a common choice for continuous forces. 

The equations of motion (5) were integrated using 
the velocity Vcrlet algorithm 40,41 with a time step of 
dt = 0.0005r or dt = 0.002r where r = a^/mj^ de- 
notes the time unit in Lcnnard Jones parameters. Using 
the described simulation scheme nVT, non-equilibrium 
(applied constant shear velocity in the walls) and equi- 
librium simulations were performed. 

We studied a melt of 10-bead polymer chains with 
segment number density close to Pmeit.coexC 3 = 0.61 be- 
tween two end-grafted polymer layers of different graft- 
ing densities. A temperature of ksT/e = 1.68 was used 



to reproduce the thermodynamic conditions of previous 
Monte Carlo and self-consistent field calculations. 7,37 At 
this temperature the density, p m eit,coex, of the melt in 
coexistence with its vapor of very low density is not 
too high and an efficient equilibration of the system by 
Monte Carlo or MD simulations is possible. If we re- 
duced the temperature further the density of the vapor 
would decrease further but the density of the coexisting 
melt would increase and, thus, significantly slow down 
the simulations. 




Figure 1: Sketch of the simulation setup. Segments of the free 
melt chains are shown in light grey, the brush is presented in 
dark grey, while the heads of the brush are light. The figure 
corresponds to a snapshot of the simulation box for a grafting 
density of R^ e /T, = 5.89 and wall velocity v x = la/r The 
shear direction, wall velocity and distance between walls are 
indicated. 

The grafting density S _1 denotes the number of 
grafted chains per unit area and we measure it in units of 
the end-to-end distance, R ee of melt chains in the bulk. 
R 2 ee /T, ~ 0(1) marks the cross-over from mushroom to 
brush regime. The system was simulated for R 2 e /T, = 
0.80, 1.47, 2.95, 5.89, 10.31, spanning different physical 
regimes of the grafted layer - from the mushroom regime 
for low grafting densities to dense, strongly stretched 
brushes exhibiting autophobic wetting behavior. Two 
smaller grafting density values, R 2 ee jY> = 0.27, 0.54, were 
also considered to analyze the effective boundary con- 
dition of the melt under shear. The area of the grafting 
substrate was taken to be 21cr x 18.2c, whereas we chose a 
fixed distance of D w = 30a between the substrates. The 
system was simulated typically for 1 x 10 6 MD steps for 
thermalization, while trajectories of length 2.6 to 4 x 10 6 
were used to average data subsequently. The walls were 
sheared relative to each other in the lateral direction, as 
shown in Fig. 1, where the main parameters of the sim- 
ulation are indicated. 

The shear velocity was taken in the range v x = 0— 2a /t 
(or alternatively shear rates < jr < 0.06). The equi- 
librium characteristics of the system were obtained from 
bulk MD simulations with periodic boundary conditions 



in all the spatial dimensions. From them, we extracted 
the end-to-end distance, R ee = Vv^Ie) = 3.66er, and 
self-diffusion coefficient, D = 0.05cj 2 /t. The characteris- 
tic relaxation time of chain conformations is thus given 
by r* = R 2 ee /D = 268r. 

To simulate shear we utilized the same thermodynamic 
conditions as in the equilibrium simulations but moved 
the first beads of the grafted chains with constant veloc- 
ity, v x , in the ^-direction and monitored the force, F x , on 
the grafting points (see Fig. 1). Using the natural time 
and length scales, t* and i? ee , we define a characteristic 
velocity, v* = D/R ee = 0.014<t/t. The ratio between the 
shear velocity, v x , measured in units of v* and the film 
thickness, D w , measured in units of R ee defines the Weis- 
senberg number, We = (v x /v*)(R ee /D w ) = jr* which 
measures the shear rate in units of the relaxation time of 
the chains. Our simulations span a wide range of Weis- 
senberg numbers for this polymer system: < We < 16. 
Adimcnsional units are used in the remaining sections un- 
less explicitly mentioned. Lengths are expressed in units 
of the chains' end-to-end distance, R ee , in the melt and 
velocities in terms of the Wcisscnberg number, We. 



III. RESULTS 
A. Segment density profiles 

First, we illustrate the general structure of the brush- 
melt interface in equilibrium and under shear by plot- 
ting profiles of the number density of segments belong- 
ing to grafted and free chains in Fig. 2. The graft- 
ing density decreases from top to bottom: i? 2 e /X = 
10.31, 5.89, 2.95, 1.47, 0.80. These values cover the en- 
tire regime from strongly stretched brushes to mushroom- 
like configurations. Upon increasing the grafting density 
we observe the building-up of an interface between the 
brush and the melt. For low grafting densities, corre- 
sponding to the mushroom regime of the brush, the melt 
chains reach the substrate. For intermediate grafting 
densities the brush layer starts to expel the melt chains 
and their density at the substrate becomes vanishingly 
small for R\ e jYi > 2.95. As we increase the grafting den- 
sity further, the average density inside the brush exceeds 
the liquid density, p me it, and oscillations in the density 
indicate fluid-like layering effects inside the brush. The 
layering becomes more pronounced when one increases 
the brush density and it is somewhat stronger than in 
previous simulations 32,33 because the grafting density is 
higher in our case. 

The brush chains stretch and the well-developed inter- 
face between melt and grafted chains moves farther away 
from the grafting surface and becomes narrower. For the 
highest grafting density studied, i? 2 e /£ = 10.31, there is 
only very little interdigitation between brush and melt. 
It is this interdigitation between brush and melt which 
is important for the static and dynamic properties of the 
system. 




Figure 2: Equilibrium density profiles for various grafting 
densities. From (a) to (e) the grafting densities decrease: 
Rle/T, = 10.31, 5.89, 2.95, 1.47, 0.80. Full lines indicate the 
density of segments belonging to grafted chains while the seg- 
ment density of free (melt) chains is marked by dashed lines. 
For all simulations, the melt density is close to the liquid- 
vapor coexistence density, p mc it,coc X o" 3 = 0.61 of the melt at 
the temperature fcsT/e = 1.68. 

To characterize the change of the brush-melt interface 
as a function of grafting density and shear velocity, we 
compute the overlap between grafted and free chains: 32 

R 2 f Dw 

^overlap = -T T / Pbrush(z) X p mclt (z)dz , (9) 

-< V Pmclt.cocx Jo 

where /9 Drus h and Pmeit are t ne brush and melt densities, 
respectively. The integrand in Eq. (9) is non-zero only in 
the regions of the film where there are free and grafted 
polymers and the integral quantifies the total amount 
of interdigitation. /overlap is directly proportional to the 
number of interactions between two polymer layers and 
therefore provides an insight on how the shear stress de- 
pends on normal pressure and shear rate. 32 ' 

Fig. 3(a) presents the overlap between brush and melt 
as a function of the grafting density. As we increase the 
grafting density, the overlap between melt and brush in- 
creases because there are more grafted chains at the sur- 
face. Upon increasing the grafting density further, how- 
ever, we observe that the overlap between melt and brush 
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Figure 3: Overlap integral as a function of the grafting density 
(a). Panel (b) presents the response to shear of the density 
profiles at intermediate grafting densities. The overlap inte- 
gral as a funstion of velocity is presented in (c). 



B. Orientation and structure 

Next we study conformational and orientation proper- 
ties. Three distinct phenomena determine the molecular 
orientation: The grafting substrate or the brush melt in- 
terface tend to align the extended molecules parallel to 
the substrate or interface, respectively. The crowding of 
the grafted chains at high grafting density results in chain 
stretching and a perpendicular orientation. Finally, shear 
tends to align molecules in the direction of the flow. 
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passes through a maximum and decreases. The behavior 
of the overlap between brush and melt parallels the wet- 
ting behavior, where one observes that the brush favors 
wetting of the melt at low grafting densities but leads to 
autophobic dewetting at large grafting densities. 

Panel (b) shows the response of the density profiles to 
shear for an intermediate grafting density. In qualitative 
agreement with simulations of polymer brushes in a good 
continuum or monomeric solvent (cf. Ref. [22] and refer- 
ences therein) the height of the brush decreases upon 
shear. Here, we observe the behavior for a brush un- 
der bad solvent condition immersed into a polymer melt. 
The thinning of the brush goes along with a significant 
decrease of the melt-brush interdigitation (cf. panel (c) 
of Fig. 3), i.e., a sharpening of the brush- melt interface, 
and it also correlates with the inclination of the grafted 
polymers for low grafting densities (see subsection IIIB). 
For the highest grafting density, i?ee/£ = 5.89, however, 
the density profiles are almost independent of shear ve- 
locity. The chains are so stretched that the shear has less 
impact on their conformations. 



Figure 4: Orientations of the brush-melt interface on different 
length scales, grafting densities and at fixed shear velocity, 
v x — la/r (or We=8.71). Panels (a) and (b) present the 
bond orientation of the brush and the melt of free chains, 
respectively. Part (c) shows the asphericity parameter A(Rc) 
of the radius of gyration. Full symbol correspond to the chains 
in the brush while open symbols characterize the behavior of 
melt chains. 

In Fig. 4 the orientation of the grafted layer and melt 
are shown as a function of the distance, z, from the sub- 
strate for various grafting densities. The bond orienta- 
tion parameter is defined as 



Sbond — 



3 (cos 2 



(10) 



where 9 denotes the angle between the vector connecting 
two bonded segments and the z direction. So defined, 
Sbond ~ 1 characterizes bonds that mainly orient along 
the z direction and Sbond ~ corresponds to randomly 
oriented bonds. Sbond ~ —0.5 denotes bonds with an 
average orientation perpendicular to the substrate (i.e., 
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along the shear direction, x, if we apply shear). The 
orientation on the scale of the entire molecule is described 
by the asphericity order parameter. We define it for the 
radius of gyration, Rq, and end-to-end vector, R ee , as 



Brush r 



A(R) 



Rl 



(ii) 



where Ri (i=x,y,z) denote the different Cartesian compo- 
nents of the considered vector. 

<Sbond for different grafting densities and for Weis- 
scnberg number We= 8.71 is presented in Fig. 4. Panel 
(a) presents the behavior of the grafted chains. In the 
ultimate vicinity of the grafting surface bonds align par- 
allel similar to the behavior of a melt in contact with 
a hard substrate. The behavior very close to the graft- 
ing surface is largely independent from the grafting den- 
sity. The bond orientation parameter exhibits some os- 
cillations that mirror the layering observed in the seg- 
ment density profiles and these oscillations become more 
pronounced at higher grafting densities (cf. Fig. 2). At 
very low grafting densities (i.e., in the mushroom regime) 
the bond orientation parameter rapidly decays for larger 
distances from the surface. At intermediate and large 
grafting densities, however, there is a pronounced per- 
pendicular orientation of the bond vectors in the mid- 
dle of the brush which stems from the stretching of the 
grafted chains. At the top of the brush the orientation 
parameter decays to zero. 

Part (b) of Fig. 4 shows the corresponding orientation 
of free chains in the melt. In the ultimate vicinity of the 
grafting surface we find a pronounced parallel orientation 
similar to the first segments of the brush. Otherwise, 
the magnitude of the orientation of bond vectors in the 
melt is much smaller than for brush chains. There is 
only a small parallel alignment of the bond vectors that 
arises from the shear. It is slightly less pronounced inside 
the brush than in the middle of the film and there is a 
slight increase of alignment as we increase the grafting 
density. This can be rationalized as follows: The larger 
the grafting density, the higher is the velocity gradient in 
the middle of the film and the concomitant orientation 
(cf. the discussion of the slip length in subsection HID). 
The behavior is similar to the parallel alignment at a 
polymer-polymer interface or a substrate. It signals the 
formation of a well-defined brush-melt interface at high 
grafting densities and the difficulty of the free chains to 
penetrate into the brush which leads to autophobicity. 

In panel (c) we investigate the behavior on larger 
length scales. Not surprisingly, orientation effects are 
much more developed for the radius of gyration than 
for bond vectors but the qualitative behavior is similar. 
Grafted chains with their centers of mass close to the 
grafting surface adopt pancake-like conformations that 
are more extended parallel to the wall than in the per- 
pendicular direction. The majority of grafted chains have 
their mass centers in the middle of the brush region 
around z R ee and they are extended away from the 
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Figure 5: Cartesian components of the average radius of gyra- 
tion as function of shear velocity, v x , for grafted chains (top) 
and free chains of the melt (bottom). Different lines corre- 
spond to various grafting densities as indicated in the key. 
For a Gaussian chain each Cartesian component of the radius 
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grafting surface. The effect increases for higher grafting 
densities and the further the chain's center of mass is dis- 
placed from the surface. Free chains in turn, are much 
less aligned and their parallel alignment arises from the 
shear, except for the two highest grafting densities, where 
one additionally observes a parallel alignment of chains 
at the brush-melt interface. 

In Fig. 5 we corroborate our findings by analyzing the 
behavior of the Cartesian components of the radius of 
gyration as a function of shear velocity. The data are av- 
eraged over all chains in the film independently of their 
distance from the surface. The top panel (a) shows the 
behavior of the brush. As we increase the grafting den- 
sity the chains stretch away from the grafting surface. 
At high grafting density the chain conformation in the 
brush hardly depends on shear velocity. At lower graft- 
ing density, shear reduces the stretching perpendicular to 
the substrate but elongates the chains in the shear direc- 
tion. The chain extension in the neutral direction, y, is 
largely unaffected by shear but decreases as a function of 
grafting density. This effect is particularly pronounced 
for the rather short chains, N = 10, utilized in our 
simulations because the chain extension in the strongly 
stretched state is not very much smaller than the con- 
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Figure 6: Average angle of the radius of gyration for the brush 
layer as a function of wall velocity for different grafting densi- 
ties. O — 90 corresponds to a brush completely perpendicular 
to the wall. 

tour length and, thus, stretching reduces the fluctuations 
of the chain extension perpendicular to the stretching di- 
rection. For longer chains, however, we would expect a 
decoupling between the chain extension parallel and per- 
pendicular to the substrate and, consequently, a smaller 
influence of the grafting density on R y G . 

The behavior of the free chains is presented in panel 
(b). The dependence of the chain dimensions on the 
grafting density is rather weak. At vanishing shear ve- 
locity and high grafting density the chains are slightly 
aligned parallel to the brush- melt interface (ir-y-plaiie) . 
Fig. 4 suggest that similar orientation effects would also 
arise at the substrate for low grafting densities but there 
are very few melt chains in the film whose center of mass 
is in the vicinity of the grafting surface. Upon increas- 
ing the shear rate, we observe an orientation of the free 
chains in the shear direction. 

It is of interest to study to what extend the conforma- 
tional changes of brush and melt chains arise from ori- 
entation of the chains or from deformation of the chains. 
In Fig. 6 the average angle of the end-to-end vector of 
grafted chains with respect to the grafting surface is 
shown. This is defined as: 

= 90° — arccos(-R ee • z/R ee ), 

where R ee is the mean end-to-end vector and z is the 
direction perpendicular to the wall. Without shear, the 
average x and y coordinates of the brush segments coin- 
cide with the location of the grafting point in the x-y- 
plane and this configuration corresponds to an angle of 
90°. Upon shearing the film the chains tilt and the aver- 
age angle decreases. The angle is smaller, the larger the 
shear rate and the lower the grafting density are chosen. 
For small shear rates the angle decreases roughly linear 
with the shear velocity, v x , and saturation effects become 
apparent once the angle drops below 45°. 

In Fig. 7 (a), the orientation of the free chains in the 
middle of the film is shown for different length scales. 
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Figure 7: (a) Orientation of the free chains in the center of the 
film versus Weissenberg number. The different panels present 
the bond orientation parameter Sbond (top), the asphericity 
parameter for radius of gyration Rg (center), and end-to-end 
vector R ee (bottom) for various grafting densities as indicated 
in the key, respectively, (b) Deformation of free chains in 
the center as a function of the adimensional shear rate. The 
theoretical results by Bruns 43 is indicated with dashed line. 
An example of the topological arrangement of the chains is 
shown in the inset. 



The alignment with the shear direction is the more pro- 
nounced the larger the length scale and, analog to the 
behavior of the brush, the orientation effects depend 
roughly linearly on the shear rate. 

This orientation effect contrasts with the deformation 
of the chain which can be measured by the change of the 
mean square radius of gyration. In Fig. 7 (b) we find 
that the coil deformation depends quadratically on the 
shear rate for small to intermediate shear velocities and 
a saturation appears for the highest shear rates. The 
data are averaged over the central region of the film con- 
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sidcring a width of lAR ee . The shear rate is also evalu- 
ated there to assure a linear velocity profile, for a precise 
calculation of the shear rate. As showed in Fig. 7, the 
curves for different grafting densities collapse, except for 
the S _1 i?g e = 2.95. This is expected due to the fact that, 
for this case, the melt density is not exactly per 3 = 0.61, 
as mentioned previously. The quadratic behavior finding 
is in accord with the theoretical prediction 43 : 



RUir) 



(o)(i 



{yrf 
8508 



(12) 



where -Rg(O) is the radius of gyration at shear rate. 
However, we see a difference with the slope predicted 
theoretically in the context of the Rouse model 44 . We 
attribute this to the topological interactions of the chains 
even though they are not long enough to be entangled 
(sec an example taken from the simulation in the inset of 
7(b)). 



C. Viscosity of the melt 

In Fig. 8 the velocity profiles of the melt under steady 
shear arc presented. From the velocity profile of the melt 
in the center of the film we extract the shear viscosity, 77, 
of the melt bulk phase according to 



n 



IVvJ ' 



where (F x ) denotes the mean force of the walls in the 
shear direction and A is the area of the substrate. Vu 
is the velocity gradient at the center. We use the linear 
velocity profile in the middle of the film (cf. Fig. 8) to 
extract the gradient according to | Vv x | = D v f 2 +b > wnere 
D w indicates the film thickness and b the slip length (see 
subsection HID ). 

The simulations were performed at constant volume 
and therefore the final value of the melt density in the 
middle of the film is unknown ahead of time and depends 
on how much the brush layer squeezes the melt. Conse- 
quently, we performed a first set of runs where the den- 
sity in the melt slightly deviated from the target value, 
per 3 = 0.61, and then adjusted the number of chains in 
the simulation cell. Alternatively, we could have simu- 
lated at constant normal pressure. The results of the 
simulations with constant density in the middle of the 
film are shown in figure Fig. 9 for the grafting densi- 
ties R 2 ee /Y> = 2.95, 5.89 and 10.31. We observe shear 
thinning typical of non-Newtonian fluids in the studied 
range of velocities which becomes very pronounced at 
even higher shear rates. In the regime of very low shear 
rates, the data have rather large error bars because of 
the very high fluctuations of force in the shear direction. 
The error bars shown in the figure are obtained from a 
block average analysis 40 of the wall force, F x . This low 
signal-to-noise ratio makes it difficult an extrapolation to 
the viscosity at vanishing shear rate. 




Figure 8: Velocity profiles as function of wall velocity for: 



(13) (a) grafting density Rt e /H = 5.89 and (b) i?L/£ 



0.80. 

The wall (and brush layer) velocity is indicated with dashed 
horizontal lines. Dashed-dotted lines show the linear extrapo- 
lation of the far-field velocity. In the first case stick boundary 
condition is observed for the fluid in the melt-brush interface. 
The latter case corresponds to a finite slip boundary condition 
in which the fluid never reaches the brush layer velocity. 



A more reliable estimate of the viscosity at zero shear 
rate can be obtained from the Green-Kubo formula. 45 
To this end we performed simulations of the bulk system 
with periodic boundary conditions in three dimensions at 
density Pmcit.coexC 3 = 0.61. The viscosity can be calcu- 
lated from equilibrium simulations as the time integral of 
the pressure tensor correlation function: 45 



VGK{t) 



V 



dt' (P a ' j (t r ) -P Q/3 (0)), 



(14) 



where V is the volume of the system and the pressure 
tensor P a/3 is given by: 



pap _ _ 



(15) 



Since the pressure tensor is a collective quantity 10 in- 
dependent simulation runs with 4 x 10 6 steps each (dt = 
0.002r) were needed for an accurate result. The average, 
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Figure 9: Viscosity of the melt as a function of the shear 
velocity for the studied grafting densities. The melt density, 
indicated in the legends was set to p = 0.61a~ 3 . The inset 
shows the viscosity as a function of t as given by the Green 
Kubo formula. 

VGK(t) in Eq. (14), over the 10 independent runs is dis- 
played in the inset of Fig. 9. Our final estimate for the 
bulk viscosity is rycK = 2.5(3). Within the error bars this 
value is compatible with the extrapolation of the zero- 
shear-rate viscosity from the non-equilibrium simulation. 
It seems that the estimates of the non-equilibrium sim- 
ulations are systematically lower than the results from 
the bulk fluctuation analysis. If there really is a sys- 
tematic discrepancy, it might be attributed to the rather 
small film thickness in the shear simulations. Especially 
at high grafting densities, where the deviations are the 
largest, the "bulk-like" middle of the film is rather narrow 
and the brush-melt interface tends to align the chains in 
the direction of the shear resulting in a somewhat lower 
estimate for rj. 



D. Effective boundary condition and slip length 

In order to exploit polymer brushes as soft, deformable 
surface coatings in micro- and nanofluidic devices two key 
characteristics of the brush have to be controlled: (i) the 
wetting properties of the liquid on top of the brush and 
(ii) the boundary condition which describes the velocity 
field of the fluid in the vicinity of the surface. Typi- 
cally, one encodes the microscopic structure and dynam- 
ics at the surface in a single parameter, the slip length 
6, which is defined as the distance, z = —6, behind the 
substrate where the extrapolation of the linear, far-field 
velocity profile, v(z), attains the substrate velocity, v x , as 
assumed by a macroscopic no-slip boundary condition. A 
positive value of the slip length, 6, implies that the fluid 
velocity at the solid surface docs not reach the substrate 
velocity, i.e., there is slip at the grafting surface, and the 



extrapolated velocity profile reaches the surface velocity, 
v x , inside the solid substrate. 

The velocity profiles in Fig. 8 also provide information 
about the hydrodynamic boundary condition. The ve- 
locity of the brush segments coincides with the velocity 
of the wall, v x , of course, and it is indicated by horizon- 
tal lines. For the free chains of the melt we observe two 
different behaviors as a function of grafting density: 

For medium to higher grafting densities (i? 2 e /E = 
2.95, 5.89, 10.31) the liquid flow on top of the brush is 
described by a stick boundary condition, in which the 
melt reaches the wall (brush) velocity for all shear rates 
studied. 46 For instance, in Fig. 8(a) the velocity profile 
of a liquid over a densely grafted brush, i? 2 e /£ = 5.89, is 
shown for different shear rates. As one can observe, the 
boundary velocity is, in fact, attained before the sub- 
strate and free chains that are inside the brush are com- 
pletely entrained. This corresponds to a negative slip 
length and characterizes the hydrodynamic thickness of 
the brush. 22 
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Figure 10: Relative slip velocity of the melt as a function of 
shear velocity for those grafting densities that have partial 
slip boundary conditions. 

At lower grafting densities (i? 2 e /E = 
0.27, 0.54, 0.80, 1.47), however, the boundary con- 
dition of the melt changes toward a finite slip. The melt 
never reaches the velocity of the grafting substrate, and 
the slip at the substrate increases as we decrease the 
grafting density. In Fig. 10 the slip velocity of the melt at 
the substrate is shown as a function of shear velocity for 
the smallest grafting densities. The boundary condition 
for these grafting densities that results in partial slip was 
also studied as a function of time. Often experiments 
report stick-slip boundary condition 27 , where the fluid 
motion switches from complete stick to complete slip in 
the course of time. We did not observe this phenomenon 
in the time window considered in the simulations. The 
system rather exhibits a steady partial slip as observed 
from the time dependence of the shear force. 

Figure 8(b) displays the velocity profile for various 
shear velocities. Although there is a notable microscopic 
slip at the substrate the concomitant slip length, b, ex- 
trapolated from the far-field velocity profile is very small 



and negative because the velocity profile strongly devi- 
ates from the linear behavior assumed in the extrapo- 
lation of the far-field velocity. Thus, microscopic slip 
velocity and distortion of the velocity profile at the sub- 
strate compensate each other and the effective boundary 
condition corresponds to no slip. 

In the limit of extremely low grafting densities 
(i2g e /S = 0.27), we observed a positive slip length (see 
inset of Fig. 12): Although the liquid wets the solid sub- 
strate, in the limit of no grafted chains, -R 2 e /£ = 0, there 
is complete slip at the substrate (i.e., b — > oo) because it 
is perfectly flat. 

We performed simulations of a single chain grafted to 
the wall to analyze the influence and spacial extension of 
the perturbation in the melt fluid velocity. The capacity 
of a single chain to drag is isolated in this way. The com- 
parisson of the mean fluid velocity in the shear direction 
is presented in Fig. 11 for a wall velocity of We=8.71. 
The system of coordinates was taken to be in the wall, 
therefore a smaller velocity means that the grafted chain 
is dragging the fluid along. Far from the grafted chain 
a high velocity is observed. Panel (a) shows the den- 
sity of beads belonging to the grafting chain and the de- 
cay of the fluid velocity along and also perpendicular to 
the shear direction. The asymmetry, with enhancement 
in the shear direction is due to the deformation of the 
grafted chain in this direction, as can be clearly observed 
in Fig. 11(b) for a plane close to the grafting point. 

The mean slip length, b, as a function of grafting den- 
sity and Weissenberg number is plotted in Fig. 12. b is 
rather insensitive to the shear rate, except for the lowest 
grafting density considered, i? 2 e /£ = 0.27, that shows a 
clear dependence with the shear rate. At very high graft- 
ing densities, the slip length is negative and its absolute 
magnitude is comparable but slightly smaller than the 
thickness of the dense brush (cf. Fig. 2). At high graft- 
ing densities the slip length does not strongly depend on 
the shear rate in accord with the insensitivity of the den- 
sity profiles in this regime. 

The errors in the slip length in Fig. 12 increase sig- 
nificantly for small grafting densities. Note that the 
chains are grafted randomly and for low grafting den- 
sities (mushroom regime), when the grafted chains do 
not strongly overlap, the influence of the grafted chains 
in the fluid is not homogeneous and quenched fluctua- 
tions of the grafting points become important (a limiting 
case is showed in Fig. 11 for a single grafted chain). This 
gives rise to sample-to-sample fluctuations which can be 
observed via the slight asymmetry of the velocity profile 
across the simulation cell. Therefore the error is large at 
low grafting densities at low grafting densities. 

As we reduce the grafting density the slip length in- 
creases. In this regime, b characterizes the rather in- 
tricate interplay between different factors: (i) As the 
grafting density decreases the physical thickness of the 
brush decreases leading to less negative values of b. (ii) 
Additionally, the melt chains penetrate deeper into the 
brush (cf. Fig. 3) resulting in a hydrodynamic thickness 
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Figure 11: (a) Average melt velocity in the shear direction 
(plain line) and density of beads belonging to a grafted chain 
(dashed line) . The left panel shows the profile along the shear 
direction and the right one in the perpendicular one. This sim- 
ulation was performed with one chain grafted in the surface. 
Panel (b) shows the two dimensional profile of melt velocity 
(up) and bead density (down) in a plane close to the grafting 
position 2 ~ 1.2cr. The center correspond to the position in 
which the chain is grafted. The shape and extension of the 
influence of the grafted chain in the melt velocities can be 
clearly observed. 



that is smaller than the physical thickness. 22 (iii) At low 
enough grafting densities there is microscopic slip corre- 
sponding to an incomplete entrainment of the melt by 
the brush. Eventually, the slip length will diverge as the 
grafting density vanishes because of our choice of a per- 
fectly smooth substrate that exhibits complete slip. 

These results suggest that the boundary condition of 
the liquid flow on top of the brush can be controlled by 
changing only the grafting density, which is readily acces- 
sible in experiments. This offers opportunities for study- 
ing flows whose boundary conditions are purely entropi- 
cally tunable without any change in the chemical nature 
of the macromoleculcs or the substrate. 
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Figure 12: Slip length as a function of grafting density for 
various grafting densities, that give rise to stick and partial 
slip boundary conditions. The points indicated with open 
circles correspond to runs made only for We= 8.71 and the 
dashed line is a guide for the eyes. The slip length 6 versus 
Weissenberg number is presented in the inset. 

IV. CONCLUSIONS 

In this work, the interface between a brush and a melt 
was investigated in equilibrium and under shear as a func- 
tion of the sliding velocities of the walls and for a wide 
range of grafting densities of the brush layer. The poly- 
mers of both, brush and melt, were taken to be physically 
identical and described by a coarse-grained bead-spring 
model. Polymer segments attracted each other corre- 
sponding to a bad solvent below the O-temperature. We 
used the DPD thermostat to duly account for the hydro- 
dynamic interactions in the system. 

The building-up of the brush-melt interface was ob- 
served and characterized upon increasing grafting den- 
sities: As expected, free melt chains were progressively 
expelled from the polymer brush layer as we increase the 
grafting density. The melt chains orient parallel to the 
brush-melt interface and bond vectors align much less 
than the end-to-end vector. This behavior is similar to 
what is observed at interfaces between immiscible poly- 
mers but here it occurs at an interface between chemi- 
cally identical species and is solely due to entropic effects. 
The pronounced structural changes at the brush-melt in- 
terface also suggest that there is a free energy cost associ- 



ated with this interface which will lead to autophobicity, 
i.e. the dewetting of the melt of identical chains from the 
brush. This was confirmed by setting up a droplet on the 
top of the brush layer which proved to be stable exactly 
for the same thermodynamic parameters than those used 
in the shear simulations. The behavior of droplets on 
this autophobic substrate will be the subject of a future 
publication. 

We complement the equilibrium information with the 
response of the brush-melt interface to shear. From the 
gradient of the velocity at the center of the film and via 
correlations of the pressure tensor in equilibrium we have 
estimated the shear viscosity of the bulk phase. The den- 
sity profiles also provide information about the hydrody- 
namic boundary condition. For low grafting densities a 
partial-slip boundary condition was found and the inter- 
face structure is strongly affected by shear. At larger 
grafting densities, however, the simulation data are well 
describable by a stick boundary condition and the brush 
is only very little affected by shear 46 . Stick-slip phenom- 
ena could not be detected during the time window of our 
simulation. Thus, a dense brush is efficient in dragging 
along the polymer melt even though thcrmodynamically 
it would dewet (at coexistence pressure). This is a first 
step towards describing the motion of droplets on this 
soft, deformable brush substrate. One key result of our 
study is the ability to change the thermodynamic prop- 
erties and hydrodynamic boundary conditions of a melt 
in contact with a substrate by changing only the grafting 
density which is an experimentally accessible parameter. 
We note that the cross-over towards autophobic behavior 
and from slip to stick boundary conditions occurs when 
the chains of the brush begin to overlap, i.e. at rather 
moderate grafting densities. 
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